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ABSTRACT 

Magnetic connections to the plunging region can exert stresses on the inner edge of an accretion disk 
around a black hole. We recompute the relativistic corrections to the thin-disk dynamics equations 
when these stresses take the form of a time-steady torque on the inner edge of the disk. The additional 
dissipation associated with these stresses is concentrated relatively close outside the marginally stable 
orbit, scaling as r~^/^ at large radius. As a result of these additional stresses: spin-up of the central black 
hole is retarded; the maximum spin-equilibrium accretion efficiency is 36%, and occurs at a/M = 0.94; 
the disk spectrum is extended toward higher frequencies; line profiles (such as Fe Ka) are broadened 
if the line emissivity scales with local flux; limb-brightening, especially at the higher frequencies, is 
enhanced; and the returning radiation fraction is substantially increased, up to 58%. This last effect 
creates possible explanations for both synchronized continuum fluctuations in AGN, and polarization 
rises shortward of the Lyman edge in quasars. We show that no matter what additional stresses occur, 
when a/M < 0.36, the second law of black hole dynamics sets an absolute upper bound on the accretion 
efficiency. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — line profiles — 
polarization — relativity 
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1. INTRODUCTION 

Early work on black hole accretion disks pointed out 
the possibility that magnetic stresses might exert a torque 
on the inner parts of the accretion disk (Page & Thorne 
1974, Thorne 1974, Ruffini & Wilson 1975, King & La- 
sota 1977). However, in virtually every recent account of 
the dynamics of accretion disks around black holes, it has 
been assumed that there is no stress at the disk's inner 
edge, which should occur very close to the radius of the 
marginally stable orbit, rms- That this should be so was 
variously argued on the basis that the plunging matter in 
the region of unstable orbits has too little inertia to affect 
the disk, or rapidly becomes causally disconnected from 
the disk, or that such stresses were due to relatively weak 
transport processes that could not compete with the large 
gravitational forces pulling matter away from the disk. Re- 
cently this view has been questioned (Krolik 1999) on the 
basis that magnetic fields are the likely agent of torque in 
accretion disks (Balbus & Hawley 1998). If this is so, and 
their strength in the plunging region is what would be ex- 
pected on the basis of flux-freezing, they should be strong 
enough in that zone to both make the Alfven speed rel- 
ativistic (postponing the point of causal decoupling) and 
exert forces competitive with gravity. 

If matter inside the marginally stable orbit does, indeed, 
remain magnetically connected to the disk, it can exert a 
sizable torque on the the portion of the disk containing 
the field-line footpoints. Gammie (1999) has shown that, 
within the confines of a highly-idealized model of inflow dy- 
namics, this torque can considerably enhance the amount 



of energy released in the disk. 

In fact, even if there were no continuing accretion, field 
lines attached to the event horizon of a spinning black hole 
and running through the disk could exert torques of a very 
similar character (Blandford 1998, D.M. Eardley, private 
communication). We will call this situation the "infinite 
efficiency limit." 

A corollary of torque on the inner edge of the disk is 
an increase in the outward angular momentum flux. In a 
time-steady state, this additional angular momentum flux 
must be conveyed by additional stress. Additional local 
dissipation must accompany the additional stress. It is 
the principal object of this paper to compute how this dis- 
sipation is distributed through the disk, and examine the 
consequences for observable properties. 

Time-steady torques at r,„s are not the only way that 
energy may be transmitted from the plunging region to 
the disk — the torque may be variable, it may be delivered 
over a range of radii, and there may be radial forces ex- 
erted that carry no angular momentum. However, in this 
paper, we will restrict our attention to this simplest pos- 
sible case. 

2. THE RELATIVISTIC CORRECTION FACTORS 

2.1. Dissipation as a function of radius 

Novikov & Thorne (1973) and Page & Thorne 
(1974) showed how the surface brightness and vertically- 
integrated stress in the fluid frame for a time-steady, geo- 
metrically thin, relativistic accretion disk could be written 
as the Newtonian forms multiplied by correction factors 
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that approach unity at large radius. In the notation of 
Page & Thorne (1974), conservation of angular momen- 
tum is given by 



dr\ ^ BdVt/dr' 



(1) 



where r is the Boyer-Lindquist radial coordinate, is the 
conserved specific angular momentum of a circular orbit 
at radius r, / is a function of radius defined such that the 
flux at the disk surface in the fluid frame F = Mof / (47rr), 
and Mo is the rest-mass accretion rate. As usual, $7 is the 
angular frequency of a circular orbit at radius r. We also 
follow Novikov & Thorne (1973) by defining four auxiliary 
functions: 



B{x) = l + a,/a;^/2 
C{x) = l-3/a; + 2a*/a;^/^ 
D{x) = l-2/x + al/x^ 
F{x) = l-2a^/x^/'^ + al/x'^, 



(2) 
(3) 
(4) 
(5) 



with X the radius in units of Tg = 



GM/c^ and a* the 
dimensionless black hole spin parameter. In the usual ap- 
proach, the boundary condition on / at the radius Vms of 
the marginally stable orbit is fms = 0. The appropriate 
boundary condition when there is non-zero stress at r^s 
is 

f - 3 Ae 

~ 2 ^ r^/^ ' ^ ' 

where Cms = C'(^ms)) and Ae is the additional radiative 
efficiency relative to the one computed in terms of the 
binding energy at Vms, eo, so that e = Ae-t-eo- This choice 
of fms ensures that the integrated additional dissipation 
matches Ae, and corresponds to a stress 



W;{rms) = 



AeMo 



(7) 



We refer to a disk with Ae = as a "Novikov-Thorne 
disk." 

Using this boundary condition, the locally generated 
surface flux becomes 



F{x) 



3 GMMq 
Sit 



3/2^1/2. 



C{x)x 



(8) 



where R^'^{x) is the expression found by Novikov & 
Thorne (1973). The standard relativistic correction fac- 
tor i?^"^ goes to zero as x approaches Xms from above, 
so that F (when the inner-edge stress is zero) peaks well 
outside the marginally stable orbit. By contrast, the ad- 
ditional dissipation due to a torque on the inner edge is 
concentrated very close to rms, and is non-zero at the in- 
ner edge. The degree of concentration can be quantified 
by measuring ri/2, the radius within which fifty percent 
of the radiation is emitted: the half-light radius is a factor 
of a few smaller for torque- driven flux than for Novikov- 
Thorne flux. Figure 1 shows the half-light radius for a 
Novikov-Thorne disk, rjyj; and an infinite-efflciency disk, 
2 ) as a function of a* . 



In the limit of infinite efficiency or zero accretion rate, 
MqC remains finite, so the first term in equation [8] domi- 
nates; in this case, the flux scales as r"'^/^ at large r rather 
than as as in the standard thin disk. The expression 
for the surface flux becomes: 



F(~)(x) 



3 



L 



3/2^1/2 



2 rgKTx"^/^ Lsdd C{x) 



(9) 



where LEdd is the Eddington luminosity and Kt is the 
Thomson opacity per unit mass. 
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Fig. 1. Half-light radii for the Novikov-Thorne disk 
(dashed line) and infinite-efliciency disk (solid line). Also 
plotted are rms (dotted line) and the region inside the 
horizon (shaded). 

The angular momentum conservation equation corre- 
sponding to equation [8] is 



/ 



dzTr4,{z) 



Mp^Kix) 
277 



3/2^1/2. 
Xms ^ms 

D{x)x^/^ 



+ RP{x) 



(10) 

where R^^ is the torque correction factor (Novikov & 
Thorne 1973, Page & Thorne 1974) in the notation of Kro- 
lik (1999). 

To clarify the meaning of the extra dissipation, wc will 
write down equation [8] pretending that gravity is purely 
Newtonian: 



F(^)(r) 



3 GMMq 




(11) 



where rj„ is the disk inner edge. The first term in the 

bracket is the usual Shakura-Sunyaev (1973) correction 
factor, while the second term is derived from the extra 
torque at the inner edge. This equation (actually first de- 
rived by Popham & Narayan 1993) never applies in the rel- 
ativistic case, but can apply, for example, to a disk around 
a star where a torque is exerted by the spinning magne- 
tosphere or through a boundary layer, or to a thin disk 
surrounding a different disk solution, such as an ADAF, 
where a torque is exerted by the flow inside the transition 
point. 



2.2. Returning radiation 

To find the surface brightness distribution of the disk as 
seen by distant observers, it is necessary first to correct the 
intrinsic surface brightness due to local dissipation for the 
additional energy supplied by photons originally emitted 
at a different radius, but returned to the disk by gravity. 
In the conventional picture, this is a small correction (Cun- 
ningham 1976). Here, however, because so much more of 
the energy is released deep in the relativistic potential, it 
can be a much greater effect. 

To compute the additional returning radiation, we fol- 
lowed the method developed by Cunningham (1976), with 
a few modifications. The numerical method is described 
in Agol (1997). We compute the flux transfer function, 
Tj, by following photons emitted from each radius that 
return to the accretion disk, assuming the disk surface 
is flat and the radiation is isotropic in the fluid frame. 
We ignore the stress carried by these photons (i.e., we set 
Tg = in Cunningham's parlance). We also assume that 
any radiation that returns to the disk inside r^s is cap- 
tured by the black hole — this radiation will be advected or 
scattered inwards by the inflowing gas, which has a large 
inward radial velocity. Finally, we (temporarily) assume 
(as does Cunningham) that the radiation returning to the 
disk is absorbed and thermalized before being reemitted; 
this assumption is probably not appropriate in practice, 
but greatly simplifies computation of the transfer function 
since in this approximation Tf is independent of frequency. 
We will discuss later how breaking this assumption may 
change the spectrum. 

In Figure 2a, we plot versus radius the fraction of emit- 
ted radiation which returns to the disk outside rms^ which 
enters the black hole or returns to the disk inside r^s , and 
which reaches infinity directly, for the cases a, = 0.9999 
and a, = 0. The fraction reaching infinity and returning 
to the disk are nearly independent of the black hole spin. 
The fraction returning to the disk is greater than 10% 
for r < 6rg, so when the emitted energy is concentrated 
inside this radius (see figure 1), then returning radiation 
will play an important role in modifying disk spectra. For 
r <j l-5rg, less than half of the radiation reaches infin- 
ity directly - most returns to the disk. For a, — 0, the 
fraction of radiation which is captured by the black hole 
or returns inside r„is increases since r^s is so large; this 
fraction never exceeds 8%. 

The fraction of returning radiation integrated over all 
radii, /rot (as measured at infinity), is shown (as a func- 
tion of a*) in Figure 2b (dashed lines) for two limiting 
efficiencies: e — cq (/"ot) and e = oo {f^t). The fraction 
for any other efficiency can be found by taking a linear 
combination of the fraction for these two efficiencies 

Act = (/°teo + /,SAe)A. (12) 

As can be seen from the figure, /ret is relatively small for 
a* = 0, even for e = oo. However, /ret grows quickly with 
increasing a*. The primary reason for this is that rms 
shrinks with increasing a*, so that relativistic effects on 
the photon trajectories become more important. Trajec- 
tory curvature is especially strong for those photons com- 
ing from small radii whose initial direction would carry 
them over the black hole. When a* and Ae are compar- 
atively large, up to 58% of the energy due to the extra 
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dissipation ends up striking the disk. 
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Fig. 2. (a) Fraction of locally emitted fiux which reaches 
infinity (solid lines), returns to the disk (dotted lines), and 
enters the black hole or returns inside r^s (dashed lines) 
versus radius. The heavy lines are for a* = 0.9999 and the 
lighter lines for a* = 0. (b) Fraction of energy that returns 
to accretion disk integrated over radius (dashed lines) or 
that is absorbed by the black hole (solid lines) as a function 
of spin. The heavy lines are for an infinite efficiency disk, 
while the lighter lines are for a Novikov-Thorne disk. 

The enhanced dissipation near r.ms also leads to an in- 
crease in the fraction of captured photons. We have com- 
puted this fraction integrated over all radii, fsH, using 
the same general relativistic transfer code just described. 
Our results for this effect are also illustrated in Figure 2b, 
again for e = eo,oo. The e = ep results agree with Thorne 
(1974). Equation applies to fsH as well. The fraction 
of locally generated radiation that ultimately escapes from 
the disk to infinity is simply fesc = 1 — fsH- Radiation 
that returns to the accretion disk we assume is reradiated 
isotropically and locally, and thus eventually reaches infin- 
ity or the black hole. We fold these multiply reprocessed 
photons into the final result. The nominal accretion ef- 
ficiency, e, is then multiplied by fesc to find the actual 
radiative efficiency of the ffow. The largest Jbh is 0.15, 
achieved for e — > cxd and a, — > 1. 

The black hole bends the radiation back to the disk so 
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that an observer on the disk sees the far side of the disk 
as a mirage above the black hole, which peaks in bright- 
ness within a few Tg of the disk plane. The flux at large 
radius then scales as H/r~^, where H represents the flux- 
weighted height of the image above the disk plane. The 
ratio of the returning radiation to locally generated ra- 
diation, Rret{f'T varies as a function of radius. For 
a Novikov-Thorne disk, Rret is infinite at r^s, then de- 
creases rapidly, asymptoting to a constant for r ^ lOr^. 
In the case of an infinite efficiency disk (with Mq = 0), 
the locally generated surface brightness scales as r~^/^ at 
large radius, while the returning radiation scales as r~^, 
so Rret diverges as r^/^ at large radius. For finite Ae, the 
returning flux may dominate at intermediate radii; how- 
ever, at large radius, Rret asymptotes to a constant due 
to the fact that for large enough radius both returning 
and locally generated flux scale as r"^. For r ^ lOvms 
and e < 1, Rret differs by at most 25% from the value 
at r = oo. We computed Rret{£,0'*,oo) as a function of 
a* and e; this function is shown in Figure 3. Fitting for- 
mulae for this quantity are given in the appendix; these 
formulae can be used to compute the returning flux at 
large radius for arbitrary a,,e. As can be seen in Figure 
3, the returning radiation can be a significant fraction of 
the locally generated radiation, and may therefore be im- 
portant for construction of disk atmospheres. Returning 
radiation can also lead to significant fluctuations on the 
light- crossing time, as discussed in §3.3. 



However, when there are stresses at Vms, angular mo- 
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Fig. 3. Contour plot of Rret at large radius. Solid con- 
tours are shown with spacing of 0.2, from Rret = at the 
bottom to Rret = 1.8 at the top. The dotted line is eo; the 
dashed line is the contour of Rret = 1- 



mentum is transferred from the matter inside 



to the 



disk. This reduces the accreted angular momentum by an 
amount Cms^ 



"I'tis^ where 



n /^l/2p— 1a, 



''ms-^ms'~-rns ^ ms 



(13) 



when all the energy liberated in the plunging region is de- 
livered to the disk in the form of work done by torque. 
^ms = S^/2Ae when a* = 0, falling towards a/SAc when 
o* approaches one. Thus, the rate at which black holes are 
spun up is substantially reduced relative to what would be 
expected in the conventional picture. Surprisingly, even 
when the black hole is initially spinless, it can be spun 
backwards when e > 1 — l/\/2 0.29! 

Considerations of black hole spin- up also place an upper 
boimd on the possible increase in efficiency due to torques 
on the disk. By the second law of black hole dynamics, the 
area. A, of the black hole must increase with time; that is. 



dA _ dA dM dA dJ 

'dt ~ dM~dr 'dj'dt ^ ^' 



(14) 



Since dM/dt and dJ/dt both depend on a, and e, this con- 
straint can be changed into a constraint on e as a function 
of 0*. For a* < 0.3584, there is a maximum achievable 
efficiency 



= 1- 



^1/2 3/2 



I + a.Ml - 2(1 + yr^) 



(15) 



Note that ^max = 1 for a* = 0, because, of course, there is 
no spin energy to tap. For a* ~ 0.3584, the denominator 
equals zero, so tmax diverges; above this critical spin, the 
decrease in angular momentum dominates the change in 
surface area, eliminating any upper bound on emax- When 
accreted radiation is included, trnax increases slightly. We 
plot tmax in Figure 4. 
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3. CONSEQUENCES 

3.1. Black hole growth and spin-up (or spin-down) 

Accreting matter enters the black hole with a certain 
amount of angular momentum, changing the spin of the 
black hole. When there is no stress at the marginally sta- 
ble orbit, the angular momentum absorbed per unit rest 
mass accreted is exactly the specific angular momentum of 

the marginally stable orbit, Lj^^ = MFmsCmV^Xms (here 
we use conventional relativistic units in which G = c = 1). 
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a. 

Fig. 4. Plot of the maximum achievable efficiency (emax) 
vs. spin due to the limit imposed by the second law of 
black hole dynamics with (solid) and without (dotted) the 
effects of radiation. 

As Thorne (1974) showed, when a* approaches one, the 
angular momentum of the black hole is also affected by 
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photon capture. Most of the photons emitted close to rms 
directed against the sense of black hole rotation are cap- 
tured by the black hole, whereas fewer of the prograde 
photons fall into the hole. This extra negative angular 
momentum prevents it from spinning up all the way to 
a* = 1. Under the assumption of isotropic radiation in 
the; fluid frame (and, of course, zero stress at rms), Thorne 
estimated that the maximum achievable a* ~ 0.998. To 
describe this effect in our context, we again normalize to 
Ll^g, so that the photon "reverse torque" per unit accreted 

mass is £-y = —{Ll^gMo)~^{dJ/dt)rad, where the notation 
is adapated from Thorne (1974). 

Combining the effects of mechanical torque and photon 
capture, we find that the net rate of change of the black 
hole's angular momentum is 



^ = LlMl - Cms 



(16) 



Because both and Cms depend on the state of magnetic 
coupling, as well as on a*, it is no longer possible to speak 
of a definite upper bound on the attainable black hole spin. 
Rather, one can instead define the accretion efficiency, 
^eq{o.*), at which da,:/dM ~ for a given a*; for e > e^q, 
the black hole is spun down due to accretion. To compute 
Ceq, we write -L^ = {dJ / dt)rad/ {MoL\^s) = J[ + eJ{ and 
[dM/dt)rad/Mo = M[ + eM^ (we give fitting formulae for 
these functions in the appendix). Then, the equilibrium 
efficiency is given by: 

2a41 + M[) - u^r'n^B,,,,^ - Ljjl + J[) 
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Fig. 5. Plot of Ceq with (solid) and without (dashed) 
effects of returning radiation. The dotted line is eo, the 
Novikov & Thorne efficiency. 

Figure 5 shows teqia*) when only magnetic and matter 
torques are included (dashed line), and when magnetic, 
matter, and radiation torques are included (solid line), as 
well as eo(a*), the efficiency of accretion (not in equilib- 
rium) when magnetic torques are ignored (dotted line). 
When accreted radiation is ignored, equation [17] simpli- 
fies to: 



= 1 - 



This limit is accurate for a* < 0.5, and only creates a sig- 
nificant error for a» > 0.9, for the radiation torque is unim- 
portant when the spin is relatively small. Some interesting 
limiting values are eeg(O) = 1 — l/\/2, and, when radia- 
tion effects are ignored, €eq{^) — 1 — which equals 
eo(l). However, when the radiation torque is included, 
Ceg is significantly reduced for a* > 0.9, and egg < eo for 
a* > 0.998, the same maximum spin found by Thorne 
(1974). The maximum equilibrium efficiency is 0.36, and 
occurs at a, = 0.94. 

A variety of spin histories is possible in this picture, for 
the efficiency is controlled jointly by the strength of the 
magnetic torques and the black hole spin. If the torque 
on the disk is always positive, the region below the curve 
Co (a,) in figure 5 is unreachable. In that case, a* = 0.998 
would still be the maximum spin achievable by accretion, 
although other spin-up mechanisms, such as black hole 
mergers or non-magnetic accretion, might permit this limit 
to be exceeded. 
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Fig. 6. Plot of the effective temperature vs. radius for 
a black hole with a» = 0.998. The temperature is normal- 
ized to To = [MoC^/(r^(jB)]^/^, where as is the Stefan- 

Boltzmann constant. The lower three curves are for a 
Novikov- Thorne disk, while the upper two are for an e = 1 
disk. The solid lines are without returning radiation, while 
the dotted lines include returning radiation. The dashed 
line shows the result of Cunningham (1976). 



3.2. Emitted spectrum 

The effective temperature is determined by the sum of 
the locally generated and returning fiux. Figure 6 illus- 
trates the effects discussed in §2, showing both how the 
intrinsic dissipation varies as a function of radius when 
there is a torque on the disk inner edge, and the total 
surface flux if one assumes that any incident radiation is 
absorbed. By comparing the curves for Ae = with the 
other curves, it is clear that the additional stress has two 
effects: the additional intrinsic dissipation creates a region 
at small radius where the effective temperature is rather 
higher than the disk could achieve otherwise; and return- 
ing radiation elevates the effective temperature at all radii, 
especially when a.^ is near unity and Ae ^ 1 or more. 

Ideally, detailed atmosphere calculations should be per- 
formed in order to ascertain the predicted disk spectrum. 
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Fig. 7. Comparison of the spectra as a function of in- 
clination angle for e = eo (dashed lines) and e = oo 
(solid lines). The other parameters are a* = 0.998, 
and Tout = 500rg. The heavy lines are for fi = 0.01 



while the lighter lines are for ^ 



0.99. 

,1/4 



The frequency 

is scaled to uq = (k/h) [L/{rgaB)]^^^ ■ The quantity 
= AnD"^ F,y{fi), where F^{fi) is the flux seen by a Eu- 
clidean observer at distance D and angle n = cos i relative 
to the accretion disk. 

with the downgoing flux of returning radiation included in 
the upper boundary condition. Interesting effects might 
well be expected due to comparable amounts of heat ar- 
riving from above as from below. Pending the completion 
of that work, we make the much simpler assumption that 
the intensity at the surface of the disk is a blackbody at the 
local effective temperature, and isotropic in the outward 
half-sphere. With that assumption. Figures 7 and 8 show 
the predicted integrated spectrum for a variety of values 
of a*, Ae, and inclination (parameterized hy n = cosi). 
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Fig. 8. Comparison of the spectra as a function of in- 
clination angle for o* = (solid lines) and a* = 0.998 
(dashed lines) with e = 1 and rout = 500rg. The heavier 
lines are for /i = 0.01, while the lighter lines for fj, = 0.99. 
Units are the same as in Figure 7. 



Figure 7 shows that for fixed luminosity and large spin, 
the efficiency of accretion can change the observed flux by 



factors of a few at different inclination angles. The angle 
dependence of the flux depends strongly on frequency 
the highest frequency radiation is concentrated towards 
the disk plane, while the lowest frequencies are radiated 
as cosi. Figure 8 shows the dependence of the spectrum 
on black hole spin for fixed luminosity and efficiency. The 
relativistic effects are much stronger for the higher spin, 
hardening the edge-on spectrum and causing strong limb- 
brightening at the highest frequencies. In contrast, the 
disk around the Schwarzschild hole is limb-darkened at 
most frequencies, and, when face-on, is brighter by a fac- 
tor of a few at the mid-range frequencies than the extreme 
Kerr hole. These effects may also impact the profiles of 
Fe Ka emission lines. If their cmissivity is proportional to 
the local flux, the enhanced flux in the inner rings of the 
disk strengthens the red wings of the lines when viewed 
more or less face-on. We plot the profiles of Ka lines for 
disks with Ae = and e = €,,q = 0.293 for a* = and 
i = 30° in Figure 9. Disks with higher spin have a smaller 
change in the shape of the iron line as a function of e be- 
cause the returning radiation is much stronger and creates 
an emissivity profile very similar to the Novikov-Thorne 
profile. Magnetized accretion may also lead to enhanced 
coronal activity immediately above the plunging region 
(Krohk 1999); if so, this would provide a physical real- 
ization for models like those of Reynolds & Begelman 
(1997), which call for a source of hard X-rays on the sys- 
tem axis a few gravitational radii above the disk plane. 
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Fig. 9. Profiles of Fe Ka lines for a* = 0, i = 30°, 
and e = = 0.293 (solid curve), e = eq (dashed curve). 
Frequency is normalized to unshifted line frequency, and 
line amplitude is normalized to the line maximum. 



3.3. Coordinated variations 

When Ae is comparable to the ordinary cSicdency, the 
inner rings of the disk radiate an amount of energy compa- 
rable to that radiated by all the rest of the disk. When, in 
addition, a* is large enough that /rot is significant, much 
of the light produced even at larger radii is reprocessed 
energy from the additional dissipation. If there are varia- 
tions in that dissipation rate, they will be reproduced — at 
appropriate delays — in the reprocessed light. A predic- 
tion of this picture is therefore that fluctuations at a wide 
range of frequencies v should all be describable as driven 
by a single source. When the fluctuations in the returning 
flux are small compared to the mean local flux (combin- 
ing both the intrinsic and the mean returning flux), the 
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relation between input and output may be written as the 
linear convolution 



dT^^{T)6L,{t-T), 



(19) 



where 5Lc{t) is the history of fluctuations in the intrinsic 
output near r^s and "ii, is a frequency-specific "response 
function" ^ that describes the distribution of relevant light- 
travel times. Note, however, that if there is a corona at 
small radii that receives a significant fraction of the to- 
tal dissipation (indeed, such a corona might receive much 
of the extra accretion energy: Krolik 1999), it will also 
drive fluctuations in the output of the outer disk in very 
much the same manner, and with a substantially identical 
response function. 

The response function ^v{t) is also predicted by this 
model. To compute this function we make several simpli- 
fying approximations: that all the returning radiation is 
absorbed; that it is reradiated in a spectrum that is locally 
blackbody and isotropic in the outer half-sphere; and that 
the radii of interest are far enough out in the disk that 
relativistic effects may be ignored. Then 



2cL, 



dr 



r3/4 



gr3/4 _|_ g-r3/4 



X I- ijl'^ -{cr/r-iy 



-1/2 



(20) 



where radius r and cr are measured in units of r*, the 
radius at which hu = kT when the flux takes its mean 

value, Tin is the innermost radius at which the return- 

-3 



mg 
r2 



flux IS (X r 
ct/(1 - ^/T 



ri 



max[rin,cT/{l - y/l - iJ? 



li^), and fji is the cosine of the in- 
clination angle. The characteristic radius r* is given by 



1/3 



(21) 



where k is the Boltzmann constant, cr is the Stcfan- 
Boltzmann constant, and is the mean value of the lumi- 
nosity emitted by the portion of the disk whose emissivity 
is oc r~^. 

Some sample response functions are illustrated in Fig- 
ure 10. All of the curves have significant tails extending 
out to ~ 10r*/c, and, almost independent of inclination 
angle, the "half-response time" (in the sense of / dr^v) 
occurs at T ~ 3r,/c. However, the peak in the response 
function becomes sharper and moves to smaller multiples 
of r*/c as the inclination angle increases. Two effects ac- 
count for this behavior. The tails are due to the fact that 
the temperature declines only as r"^/^, so that the Wien 
cut-off sets in relatively slowly. The sharp peaks at small 
lag exhibited by disks with larger inclination angle are due 
to the significant amount of disk surface that lies close to 
the line of sight for those viewing angles. 




Fig. 10. Plot of continmmi response func- 
tion of lag for different inclinations. is in units of 
fieti"*hi'^ / {2cL^) and r is in units of r*/c. 

3.4. Polarization 

In fact, the inner rings of realistic disks, whether in 
AGN or Galactic black holes, are likely to be scattering- 
dominated, so that their albedo to the returning radiation 
will be significantly greater than zero. The scattered light 
may then be polarized. 

As discussed in §2.2, the maximum altitude H above the 
disk plane achieved by any photon that ultimately returns 
to the disk cannot be much greater than a few gravita- 
tional radii. If the disk flare is small (sec §3.6 for further 
discussion), the returning photons striking the disk at ra- 
dius r must then arrive from an angle ~ H/r <C 1 from 
the disk equator. When electron scattering is the domi- 
nant scattering opacity (as is nearly always the case) , only 
those photons polarized parallel to the disk normal can 
scatter to outgoing directions near the equatorial plane 
but perpendicular to the original photon direction. The 
result is that disks viewed obliquely should acquire a small 
amount of polarization parallel to the disk axis, especially 
at the high frequencies produced predominantly in the in- 
ner rings. 

To quantify this suggestion, we have computed the disk 

spectrum, treating the locally generated radiation as a 
blackbody, and assuming the returning radiation is scat- 
tered off a semi-infinite electron scattering atmosphere 
(Chandrasekhar 1960). We assume the locally generated 
disk flux has either (1) the polarization of a semi-inflnite 
electron scattering atmosphere (Chandrasekhar 1960) or 
(2) is unpolarized. The true polarization will be mod- 
ified by Faraday rotation due to magnetic fields in the 
disk's atmosphere and absorption/emission (Agol, Blaes, 
& lonescu-Zanetti 1998), but the true answer will likely 
lie between our two assumptions. Figure 11 shows the 
flux, polarization, and polarization angle computed under 



^We use the term "response function" to avoid confusion with the relativistic "transfer function" 
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these two assumptions. The spectrum is mueh broader 
than would be predicted by complete absorption, and re- 
turning radiation can cause a sharp rise in polarization 
towards the highest frequencies. There is a rotation in the 
polarization angle since the scattered returning radiation 
tends to be polarized perpendicular to the disk plane, while 
the locally generated radiation is polarized parallel to the 
disk plane. Whichever component dominates the flux at a 
given frequency determines the strength and angle of the 
polarization. We have included all relativistic effects that 
modify the final polarization angle (Laor, Netzer, & Pi- 
ran 1990, Agol 1997). The returning fraction is largest for 
photons generated in the inner region of the disk; conse- 
quently, the highest frequencies have the largest scattered 
fraction and thus the highest polarization. In addition, 
the inner parts of the disk are strongly blueshifted, and 
the returning radiation is (weakly) Compton up-scattered 
by the bulk motion of the disk. 

0.0^ — . — . — . — , — . — . — . — . — , — . — . — . — . — , — . — . — . — . — , — . — . — . — ^ 



-0.5 - 




-1.5 -1.0 -0.5 0.0 0.5 1.0 



Fig. 11. Flux, polarization, and polarization angle as a 
function of frequency for a disk viewed with an inclination 
of = cosi = 0.2, with a = 0.998, Ae = 1. The dashed 
curves are for no returning radiation; dashed-dot for re- 
turning radiation, but unpolarized locally generated flux; 
the solid curves are for returning radiation plus polarized 
locally generated flux; and the dotted line in the top panel 
shows the flux computed assuming complete absorption, 
as in §3.2. The polarization angle, 9, is zero for E parallel 
to the disk plane. The units are defined in Figure 7. 

3.5. Bolometric Limb-brightening 

For a Newtonian disk, foreshortening causes limb- 
darkening proportional to /x = cosi, where i = is a face- 
on disk. Relativistic effects cause beaming and bending of 



the radiation towards the equatorial plane, which decrease 
the limb-darkening for disks around black holes. For large 
a* and e, the relativistic effects become so strong that a 
disk can actually become limb-brightened. In Figure 12 we 
show the bolometric disk flux as a function of inclination 
angle for the cases e = oo (dashed fine), e = 1 (dotted 
fine), and e = eo (solid fine) for a, = 0.998. The limb- 
brightening is also be dependent on frequency as shown 
in §3.2; in practice, determining this quantitatively will 
require a detailed disk atmosphere model. 



3.0 -\ 
2.5 - 
2.0 - 
a: 1.5 - 




0.5 - 



0.0 [ > > > I > > > I > > > I > > > I > > > I 

0.0 0.2 0.4 0.6 0.8 1.0 

Fig. 12. Comparison of the bolometric limb-brightening 
of a disk for o* = 0.998 with e = eo (solid line), e = 1 (dot- 
ted line), and e = oo (dashed line). The inclination angle, 
jj, = cosi is edge-on for /i = and face-on for ji = 1. 

3.6. Geometrical thickness of the disk 

When the accretion rate is greater than a small fraction 
of the Eddington rate, the innermost regions of accretion 
disks are expected to be supported against the vertical 
component of gravity by radiation (Shakura & Sunyaev 
1973). In that case, the disk's vertical thickness is di- 
rectly proportional to the ratio of the local radiation flux 
to the vertical component of gravity; that is, h oc Fr^ /R^, 
where is the relativistic adjustment to the vertical grav- 
ity (Page & Thorne 1974; Abramowicz, Lanza, & Percival 
1997). At radii large enough that the relativistic effects are 
small, but not so large as to no longer be in the radiation- 
dominated regime in a Novikov-Thorne disk, h should be 
constant. In the relativistic portion of the disk, h would 
shrink oc R^'^ / Rz if the stress at its inner edge were zero; 
additional dissipation, depending on its strength, could ac- 
tually make the disk become somewhat thicker there (cf. 
equation 8). We plot some examples of h{r) in Figure 13. 
At the inner edge of the disk, the height of the disk is 
non-zero when there is a non-zero torque, so the thin-disk 
approximation is valid only if h{r„is) < 0.1r„,s. This crite- 
rion can be translated into a limit on the extra luminosity 
due to the torque at the inner edge: 

^ C 1 ^ (x C-^/'^F^ -a^G -Ka^C^/^^ 

^Edd oXms ^ ' 

(22) 

where L = ^eMoC? . This limiting luminosity is plotted 
in Figure 14. When Ae = 0, the thin disk approximation 
breaks down if hjr ^ 0.1 where h/r is maximum. This 
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limit can in turn be expressed as a limit on the luminosity 
coMqC^, which is also shown in Figure 14. The luminosity 
upper limit for the infinite efficiency disk is much smaller 
than for the Novikov-Thorne disk since /;,/ r peaks at rms , 
while for the Novikov-Thorne disk h/r peaks at larger ra- 
dius, where its magnitude is smaller [{h/r)max occurs at 
r = 24rg for a* = and r = 7rg for a, = 1]. If either 
e-oMoC^ or AeMqC^ exceeds their respective limits, then 
the thin-disk approximation breaks down. In addition, if 
h/r is small, then the approximation of a flat disk in the 
computation of the returning radiation will be appropri- 
ate. To treat the interesting cases where L ~ Lsdd will 
require a 2-D solution of the disk equations, which is be- 
yond the scope of this work. The returning radiation will 
not affect the disk height since it diffuses through the disk 
on a thermal timescale, so there is no net flux due to re- 
turning radiation (unless the disk is warped). 



Thorne disk; solid line is for infinite-efflciency disk. 



1 .4 




Fig. 13. Height of radiation pressure-supported disk 
vs. radius with e = eo (dotted line) and e = 1 (solid 
line) for a* = 0.998. The height is normalized by ho = 
3K;Mo/(87rc). 
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Fig. 14. Upper limit on the luminosity for thin-disk 
{h/r ^ 0.1) approximation to be valid at rms for a radia- 
tion pressure-supported disk. Dotted line is for Novikov- 



4. CONCLUSIONS 



We have generalized the equations for an azimuthally 

symmetric, geometrically thin, time-steady accretion disk 
around a black hole to include the effects of a torque op- 



erating at the inner boundary, taken to be at 



Con- 



stant non-zero torque at rms causes several physical con- 
sequences that change the fundamental properties of the 

accretion flow: 

1) The flux can be expressed as a sum of the usual 
Novikov & Thorne expression plus a part due to the torque 
which scales roughly as r^^/^. 

2) The accretion efficiency has a fundamental upper 
limit due to the sec;ond law of black hole dynamics for 
a* < 0.36. For larger a*, infinite efficiency is possible in 
principle. 

3) The black hole spin can reach an equilibrium for 
a* < 0.998 since the angular momentum reaching the hole 
is smaller. Radiation can also exert a significant torque on 
the black hole, which changes the value of the equilibrium- 
spin efficiency. Above an efficiency of e = 0.36, the black 
hole must always be spun down. 

4) Since the extra emissivity is peaked at the inner edge 
of the disk, if rms is small then gravity causes a large frac- 
tion of the radiation (up to 58%) to return to the accretion 
disk. The flux of returning radiation scales as r~^ at large 
radius. Up to 15% of the radiation can be captured by the 
black hole. 

5) The extra heating within the disk will increase the 

height of the disk if it is radiation pressure-supported. 
This limits severely the luminosity at which the thin-disk 
approximation is appropriate. 

6) Doppler beaming and relativistic bending are 
strongest in the inner parts of the accretion disk where 
the extra flux peaks, so that for large a* and e, the disk 
will be limb-brightened. 

These each have multiple observable consequences: 

1) The extra surface brightness changes the locally ra- 
diated spectrum. Though the local surface brightness is 
usually not directly observable, it may be possible to map 
it using several devices: eclipse mapping (Baptista et al. 

1998) , although no eclipsing black hole X-ray binaries have 
been discovered yet; quasar microlensing (Agol & Krolik 

1999) ; or reverberation mapping (Collier et al. 1999). If 
the Fe Ka emissivity is proportional to the local dissi- 
pation, these effects can strengthen the red wing of the 
line, particularly when the spin is small. This effect may 
undercut the argument that lines with strong red wings 
came from disks around black holes with higher spin (e.g. 
Dabrowski et al. 1997). 

Several authors have used the Novikov-Thorne model to 
flt the soft X-ray spectra of galactic black hole candidates. 
Their procedure was to estimate the effective radiating 
area required to emit the observed luminosity at the ob- 
served effective temperature. On the basis of these fits 
they inferred that some black holes have rather high spins 
because the effective radiating area of a Novikov-Thorne 
disk decreases with increasing spin (Zhang, Cui, & Chen 
1997). However, for fixed spin and central mass, a disk 
with large Ae has a smaller effective radiating area than a 
Novikov-Thorne disk, mimicking the effect of greater spin. 

2) The outer parts of accretion disks can be unstable to 
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warping due to irradiation from the center (Pringle 1996). 
The minimum radius for growth of small warps is propor- 
tional to e~^; if the efficiency is much higher than that 
of a standard disk, the minimum radius may be greatly 
shrunk. Limb-brightening increases the effective efficiency 
and therefore makes the linear growth more rapid; on the 
other hand, the corresponding relative decrease in inten- 
sity away from the central disk plane may weaken this 
effect in the non-linear regime. 

3) Wavelength-dependent limb-brightening (or limb- 
darkening) introduces viewing angle-dependent biases into 
any flux-limited sample. This is a particularly strong effect 
in the context of quasar surveys because the number count 
distribution is so steep. A variety of distortions could oc- 
cur in our view of what constitutes a "typical" quasar (cf. 
Krolik & Voit 1998). 

4) Returning radiation can change the conditions for 
launching a radiation-driven disk wind (e.g. Murray et 
al. 1995). The returning radiation is reradiated locally, so 
the net vertical force depends on the frequency-averaged 
opacities of the downgoing and upgoing radiation fields. 
The radial component of the radiation force will be larger 
than for a standard disk due to the higher efficiency and 
limb brightening. 

5) Returning radiation causes the various annuli to 
"communicate" on the light crossing timescale. Fluctu- 
ations of the flux at small radii will cause only slightly de- 
layed fluctuations of the flux at larger radii, which emit at 
longer wavelengths. Indeed, exactly this sort of behavior is 
commonly seen in accreting black hole systems. For exam- 
ple, campaigns monitoring AGN have consistently found 
that continuum fluctuations are very nearly simultaneous 
afl the way from ~ 1300 A to ~ 5000 A (Clavel et al. 1991; 
Korista et al. 1995; Wanders et al. 1997; Collier et al. 
1998; O'Brien et al. 1998; Cutri et al. 1985). Comparing 
the upper bounds on any inter-band delays to the radial 
scales expected on the basis of conventional disk models, 
these observations have been interpreted as requiring a 



coordinating signal group speed of at least ~ 0.1c (e.g., 
Krohk et al. 1991; Courvoisier & Clavel 1991; Collier et 
al. 1998). 

6) The scattered component of returning radiation is 

highly polarized parallel to the disk axis at high frequen- 
cies. This polarization rise may be related to the observed 
sharp rises in several quasars (Koratkar et al. 1995, Impey 
et al. 1995). If the inner regions of the disk have a strong 
Lyman continuum in emission, then this scattered emis- 
sion edge will appear as a strongly polarized, blueshifted 
emission edge in the spectrum. In an irradiated disk atmo- 
sphere there might be an additional effect: heating of the 
upper layers of the atmosphere can cause a temperature 
inversion, which changes the sense of the polarization. We 
leave all such detailed calculations to future work. 

A question left unanswered by this work is what spin 
and efficiency we expect to be achieved by black holes in 
nature. That there is an upper limit on efficiency for an 
equilibrium spin means that a black hole with e > Ceq must 
be born with original spin, or must be spun up by accre- 
tion in which magnetic torques inside rms do not play an 
important role. In the supermassive black hole case, the 
spin may result from a merger. No one has computed the 
final spin of the resulting merger of two black holes; how- 
ever, current approximate calculations indicate that the 
final spin could be quite large (Khanna et al. 1999). The 
strength of the torque at rms depends on the strength of 
the magnetic field in the accretion disk and the geomet- 
rical thickness of the flow, which in turn depend on the 
accretion rate. This dependence will be best addressed 
with numerical simulations of MHD accretion. 
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APPENDIX 

The function iJ^et can be expressed as the sum of a part due to the Novikov & Thorne accretion rate, and a part due 
to the torque at the inner edge: 

Rret{e,a*,r) = Ro{a^,r) + Roo{a^,r)Ae. (1) 
At large radius, these functions become constant, where we have fitted them with polynomials in x = logio{^ — a*): 

i?o(a,, oo) = 0.0200 - 0.0360x + 0.0279a;2 + 0.00213.t3 - 0.00153.t'' - 0.000225a;^ 
i?oo(a*> oo) = 0.594 - 0.199a; - O.liea;^ - O.lOTx^ - 0.0373a;'^ - 0.00409x^ (2) 



Table 1 

Coefficients for fits to dM/dt and dJ/dt. 



fit 


ao 


Ol 


02 


as 


04 


as 


ae 


lOHdM/dt)°^^ 


2.402 


-11.63 


10.87 


-28.47 


-17.30 


-3.651 


-0.2731 


lQ\dMldt)Zd 


261.8 


-600.9 


449.7 


554.9 


224.9 


41.28 


2.885 


l^\dJ/dtt^, 


.5182 


1.055 


-5.149 


3.707 


2.347 


.5455 


.04527 


10\dJ/dt)^^^ 


88.53 


29.56 


-29.19 


-200.6 


-128.1 


-30.00 


-2.475 
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The formulae are accurate to better 0.6% from a* = to a* = 0.9999. 

For computing the spin evolution of an accreting black hole, it is useful to know {dJ/dt)rad = {dJ/dt)'^g^^ + A€{dJ/dt)^j 
and {dM/dt)rad = {dM/dt)^^^ + Ae{dM / dt)^^^. Note that {dM/dt)^^^ = M[ + eoM!^ and {dM/dt)'^^^ = M!^, and likewise 
for J. We have fitted these with 6th order polynomials in x. The coefficients ai of the fits (X]i=o c-i^^) given in table 
^, where we have multiplied each by 10'*. 
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